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O : Abstract 

Q-ij With LIGO having achieved its design sensitivity and the LIGO S5 strain data being available, 



< 



constraints on the relic gravitational waves (RGWs) becomes realistic. The analytical spectrum 
of RGWs generated during inflation depends sensitively on the initial condition, which is 
generically described by the index f3, the running index at, and the tensor-to-scalar ratio r. 
By the LIGO S5 data of the cross-correlated two detectors, we obtain constraints on the 
£jT)! parameters (/3, a>t,r). As a main result, we have computed the theoretical signal-to noise ratio 

(SNR) of RGWs for various values of ((3,at,r), using the cross-correlation for the given pair 
J> ■ of LIGO detectors. The constraints by the indirect bound on the energy density of RGWs by 

BBN and CMB have been obtained, which turn out to be still more stringent than LIGO S5. 

G\ ' 

PACS numbers: 04.30.-w, 04.80. Nn, 98.80.Cq 

o . 

1. Introduction 

Recently, LIGO S5 has experimentally obtained so far the most stringent bound on the spectral energy 
density of the stochastic background of gravitational waves, Qq < 6.9 x 10 around ~ 100Hz [1]. 
Generated during inflation, RWGs is of cosmological origin, and has long been investigated [2 El IH [5] , 
and, in particular, its analytical spectrum has been known [6J. It depends most sensitively upon the initial 
condition, which can be generically summarized by the initial amplitude, the spectral index /3, as well as the 
running index at- In particular, small variations of /3 and at will cause substantial change of the amplitude 
in higher frequencies [7]. The value of (3 and at are usually predicted by specific inflationary models [8J with 
possible modifications by quantum field renormalization [9J. After inflation, RGWs is altered substantially 
only by a sequence of subsequent expansions, the reheating, the radiation era, the matter era, and the 
current acceleration era |6J, essentially unaffected by the cosmic matter they encounter. As a result, 
RGWs carry a unique information of the early Universe, and can probe the Universe much earlier than the 
cosmic microwave background (CMB). Such cosmic processes, as neutrino free-streaming [10, 1 lj , QCD 
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transition, and e + e~ annihilation [12J, etc, affect RGWs less substantially than small variations of /3 and 
at around the frequency range ~ 100Hz of LIGO, and can be neglected in this study. 

Spreading over a broad range of frequency, (10 -18 ~ 10 10 ) Hz, RGWs is a major target of detectors 
working at various frequencies, including LIGO [13J, Advanced LIGO [HJ, LISA [15], EXPLORER [To], 
millisecond pulsar timing |17| , and Gauss Beam [T5], etc. The curl type of CMB polarization is only 
contributed by RGWs, measurements of which also serve as detectors [19], such as WMAP [201 [211 [221 
[23j[2lj, Planck [25] and CMBpol [26]. Prior to the LIGO S5 bound [lj, often used were the bound from 
big bang nucleosynthesis (BBN) [27] and that from the CMB anisotropy spectrum [28J. These two indirect 
bounds actually constrain the energy density f Q g (f)dln f as an integration of the spectral energy density 
Q g (f). By contrast, the LIGO S5 bound is upon Q g (f) itself, and has surpassed the LIGO S4 [29] by 
more than an order of magnitude. It is now realistic to infer from this bound some constraints on the 
initial condition of RGWs in terms of (/3, at,r). In this letter, using the strain data from LIGO S5 [TJ, we 
will derive such constraints, and compute the theoretical SNR for the analytic spectrum of RGWs with 
various at,r). 

2. Analytical Spectrum of RGWs 

In a spatially flat Robertson- Walker spacetime, the analytical mode /ifc(r) of RGWs is known [6J. The 
spectrum at the present time th is given by 

h(k,T H ) = ^k 3 / 2 \h k (T H )\, (1) 
7T 

related to the characteristic amplitude [3D], h c (f) = t#)/\/2. Here the frequency / is related to the 
wavenumber k via / = kHo/lirj with 7 ~ 1.97 for S1a — 0.73 [llj. The spectral energy density [3J [7] 

^(/) = ^(/)(^) 2 , (2) 

where H = 3.24 h x 10 -18 Hz. RGWs is completely fixed, once the initial condition is given, which is 
taken at the time Tj of the horizon-crossing during the inflation, of a generic form [7, 20, 22J 

h(k,r t ) = A^o)r5(A)2+/Hi« t ln(W (3) 

where the pivot wavenumber k corresponds to a physical wavenumber 0.002 Mpc -1 , the tensor-to- 
scalar ratio r = A^(/co)/A^(A;o) is a re-para metrization of the normalization Ah(ko) with A^(fco) = 
(2.445 ± 0.096) x 10~ 9 by WMAP5+BAO+SN Mean [22J, the index /3 is related to the index of the 
power-law scale factor during inflation a(r) oc |t| 1+ ' 3 [3j [6] and (3 ~ — 2 yields a nearly scale-invariant 
spectrum, and the running index at reflects an extra bending. Observations of CMB anisotropies have 
given preliminary results on the scalar index and the scalar running index [20, 21, 22J. So far there is no 
observation of /3 and a t , and there are only some upper bound on r [22J [231 [21] ■ l n scalar inflationary 
models, /3 and a t are determined by the inflation potential and its derivatives [8J. There might be relations 
between the tensorial indices and the scalar ones. For generality, we treat (P,a t , r) as independent 
parameters. In literature the notation n t is often used, nt = 2/3 + 4. 
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Figure 1: Left: The spectrum h c (f)yF/y/2j with r = 0.55 and at = for various j3. The dot 
line (labeled by Hl/Ll goal) is the strain sensitivity hf of single interferometers achieved by LIGO 
S5, and the dash line (Oq = 6.9 x 10 -6 ) is the corresponding sensitivity from cross-correlated two 
interferometers of LIGO S5 [TJ. Right: The spectral energy density £l g (f) with at = for various j3. 
The upper limit of LIGO S5 corresponds to the dash line in the left. 



3. Constraints on the spectral indices of RGWs 

The left panel of Fig. [TJ gives the analytic spectrum h c (f)y/F /y/2j of RGWs in the frequency range 
(40,500) Hz for various /3 in the model r = 0.55 and at = 0. The irregular oscillations in the curves of 
the analytic spectra are due to the combinations of Bessel functions implicitly contained in the analytic 
solution of RGWs [6, 7J. It is seen that a small variation in j3 from —2.0 to —1.85 leads to an enhancement 
of amplitude of h c (f) by 4 orders of magnitude around ~ 100 Hz. For RGWs to be detectable by a single 
detector with a strain sensitivity hf, the condition is |30j . 



K(f) 



2/ 



F > h 



(4) 



where the angular factor F = 2/5 for one interferometer. The dot line (labeled by Hl/Ll goal) in the 
upper part of the left of FigQ]is the single-detector strain sensitivity achieved by HI and LI of LIGO S5 [lj. 
Thus, we have plotted h c (f)y/~F /y/2f to directly compare with the strain hf. The single interferometers, 
HI and LI, of the LIGO S5 put a constraint on the index: j3 < —1.85 for the model r = 0.55 and at = 0. 

However, by the cross-correlation of two interferometers, HI and LI, of the LIGO S5, the detectability 
is much improved. Approximately, in a narrow band Af of frequencies and a duration T of observation, 
the detectability condition is schematically changed to [30] 



hc(f) 



F > 



-In 



(5) 



y/2j " ' (2TA/) 1 /4' 

where ht is the strain of single detector. For T being long enough so that (2TA/) 1 / 4 3> 1, the right hand 
side of Eq.([5]) will be reduced considerably. A detailed description of quantitative treatment is given in 
Ref.|3l]. For the case of a flat spectral energy density f2 . tne effective strain of LIGO S5, plotted in the 
dash line in left of Fig(TJ is ~ 100 times lower than that from the single interferometers [I]. This upper 
limit leads to a more stringent constraint on the index: j3 < —1.88 for the same model. This is consistent 
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Figure 2: Left: h c {f)y/F / y/2f with r = 0.55 and /3 = -2.0 for various a t . Right: Q g (f) with r = 0.55 
and /3 = —2.0 for various at- 

with the current observational result of the scalar index n s ranging over (0.97 ~ 1.2) [201 EE 122], if a 
relation n s = 2/3 + 5 is adopted, as in scalar inflationary models. 

The right of Fig. [Hgives the spectral energy density ft g (f ) that corresponds to the respective spectrum 
h c (f) in the left. By the upper limit Q, = 6.9 x 10~ 6 from cross-correlated interferometers of LIGO S5, the 
resulting constraint is f3 < —1.88, the same as from the left. Except for the model j3 = —2.0 and at = 0, 
£t g (f) is generally not flat, and a larger (3 leads to a higher amplitude of £l g (f) in higher frequencies [3J |7]. 
£l g (f) behaves approximately as £l g (f) oc / 024 for the model /3 = —1.88 and at = 0. For comparison, 
the sensitivity of LISA is plotted and has a broader frequency range. 

The left of Fig. [2] plots h c {f)y/~F / y/2j for various a t in the model of r = 0.55 and /3 = -2.0. A small 
variation in at from to 0.018 enhances the amplitude of h c (f) by ~ 4 orders of magnitude around ~ 100 
Hz. The single interferometers of the LIGO S5 puts a constraint on the running index: a t < 0.018. The 
cross-correlation of two interferometers of the LIGO S5 puts a more stringent constraint: at < 0.01. So far 
the preliminary observed result of the scalar running index a s ranges over (—0.050 ~ —0.077) by WMAP 
[221 [2U [22] . If both RGWs and scalar perturbations are generated by the same inflation, one expects at 
to be nearly as small as a s for several kinds of smooth scalar potential |8J. If so, the constraint on at by 
LIGO S5 is consistent with the results by WMAP. The right of Fig. [2] gives £l g (f) that corresponds to 
those in the left. The upper limit of LIGO S5 gives the constraint a t < 0.01, same as that from the left. 
For the model (3 = -2.0 and a t = 0.01, the slope is oc f 0A5 , not flat either. 

FigureSshows that, around ~ 100Hz, the model /3 = -2.0 and a t = 0.011 and the model /3 = -1.88 
and at = yield the same height of amplitude detectable by LIGO S5. Moreover, the slopes of O s (/) in 
the two models only differ slightly. Therefore, there is a degeneracy between the indices (3 and a t . Given a 
rather narrow frequency range, (41.5, 169.25) Hz, it is unlikely for LIGO S5 to distinguish the spectra from 
these two models. Comparatively, LISA with a much broader frequency range would have consequently a 
better chance to distinguish models with different /3 and at- 

The above examinations on detectability via comparison of the spectrum h c (f) and the strain hf are 
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Figure 3: f2 ff (/) has the same height at 100Hz for the models with f} = —2.0 and at = 0.011, and 
with P = —1.88 and at = 0, respectively. 

Table 1: The SNR for RGWs with r=0.1 for the given pair of detectors of LIGO S5 
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still qualitative. According to the method developed in Ref.|31], a more quantitative description of the 
detectability is through the signal-noise ratio 



10vr 2 



df i 2 (mi(f) 11/2 



(6) 



/ 6 P1(/)P 2 (/) 

for the given pair of detectors of LIGO, where P\(f) and P^if) are the noise power spectrum of detector, 
HI and LI, respectively [I], and 7(/) is the overlap reduction function [31]. Since the data of the strain 



sensitivity hf\ = ypjj) and /iy 2 = yPAJ) have been given [I], it is straightforward to calculate 
SNR from £l g (f) for each model. For the model Q\ = 0.73 and $7 m = 0.27, we have computed the 
corresponding SNR for various indices P, and at, listed in Table 1 for r = 0.1. The duration T in Eq.© 
for LIGO S5 is from Nov. 5, 2005 to Sep. 30, 2007 [1], i.e., T = 59961600 seconds. Clearly, greater 
values of P and at yield higher SNR accordingly. For other values of r, the corresponding SNR follows 
immediately since SNR oc r. 

3. Constraints via the energy density Vt gw 

Before LIGO S5 data is available in constraining the spectrum Q g (f), often used is the energy density 
parameter 

ff upper Jf 

*V = / , (7) 

''flow J 



as an integration of £l g (f) over certain frequency range, where the cutoffs of frequencies depend on 
specific situation under consideration. For the total energy density of RGWs in the universe, one can 
take fi ow ~ 2 x 10~ 18 Hz and f uvveT — 10 10 Hz [11J. Strictly speaking, limits coming out of this 
method do not apply to the spectrum Q g (f), and are of indirect nature. Sometimes £l g {f) and Q gw were 
used undiscriminatingly in literature. But this will be valid only under the condition that the integration 
interval df/f = din/ ~ 1 and that O g (/) be nearly frequency-independent (flat), which is not the case 
for general indices j3 and at, as has been demonstrated earlier. Whenever possible, one should distinguish 
£l g (f) and Q, gw for a pertinent treatment. Currently, two observed bounds on £l gw are available. One is 
Q gw < Qbbn = 1.1 x W^ 5 (N l/ — 3) from BBN, where N v is the effective number of relativistic species 
at the time of BBN. The abundances of light-element, combined with WMAP data, give (N u — 3) < 1.4 
|27| . so Qbbn = 1-5 x 10~ 5 [25]. This bound receives contribution from frequencies down to the lower 
limit flow ~ 10~ 10 Hz, corresponding to the horizon scale at the time of BBN |31J. Another bound 
is Q gw < ClcMBh 2 = 8.4 x 10~ 6 at 95% C.L. from CMB + matter power spectrum + Lya for the 
homogeneous initial condition of RGWs [28]. For the Hubble parameter h = 0.701 [22J [23], this is 
&CMB = 1-62 x 10~ 5 , receiving contributions from frequencies down to a much lower limit fi ow ~ 10~ 15 
Hz, corresponding to the horizon scale at the decoupling for CMB. From the theoretical side, substituting 
the analytical spectrum £l g (f) as the integrand into Eq.([7]), the resulting integral Q, gw is a function of the 
indices (3 and a t , since £l g (f) intrinsically depends on j3 and at- By this way, we can derive constraints 
on /3 and at by the bounds Qbbn and Qcaib- In carrying out the integration, we take the upper limit 
of integration f upP er = 10 10 Hz [TT]. As for the lower limit, we take fi ow = 10~ 10 Hz for BBN case, and 
flow = 10 _15 Hz for CMB case, respectively. It turns out that the integral Q gw is sensitive to the value of 
flow for very small j3 and at- 

The left of FigfJ] shows the j3— dependence of Q gw for fixed at = and r = 0.55 and 0.1, and the 
right shows the at— dependence of Q, gw for fixed /3 = —2.0 and r = 0.55 and 0.1. In FigfJ] the horizontal 
dash lines are the bounds £Ibbn and &cmb< which are close to each other. The resulting constraint on 
P is P < —1.96 for r = 0.55 and at = 0, and P < —1.98 for r = 0.1 and at = 0. The resulting constraint 
on a t is a t < 0.004 for r = 0.55 and P = -2.0, and a t < 0.005 for r = 0.1 and P = -2.0. These 
constraints on p and at by BBN and CMB are more stringent than those by LIGO S5. 
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